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We study the static and the dynamic response of coherently coupled two component Bose-Einstein 
condensates due to a spin-dipole perturbation. The static dipole susceptibility is determined and 
it is shown to be a key quantity to identify the second order ferromagnetic transition occurring 
at large inter-species interaction. The dynamics, which is obtained by quenching the spin-dipole 
perturbation, is very much affected by the system being paramagnetic or ferromagnetic and by the 
correlation between the motional and the internal degrees of freedom. In the paramagnetic phase 
the gas exhibits well defined out-of-phase dipole oscillations, whose frequency can be related to the 
susceptibility of the system using a sum rule approach. In particular in the interaction SU (2) sym¬ 
metric case, i.e., all the two-body interactions are the same, the external dipole oscillation coincides 
with the internal Rabi flipping frequency. In the ferromagnetic case, where linear response theory 
in not applicable, the system show highly non linear dynamics. In particular we observe phenomena 
related to ground state selection: the gas, initially trapped in a domain wall configuration, reaches 
a final state corresponding to the magnetic ground state plus small density ripples. Interestingly 
the time during which the gas is unable to escape from its initial configuration is found to be 
proportional to the square root of the wall surface tension. 

PACS numbers: 


I. INTRODUCTION 

Ultra-cold gases allow the realisations of multi- 
component Bose-Einstein condensates (BECs). The lat¬ 
ter are novel systems, whose behaviour is very different 
with respect to the one of single component BEC. In 
particular they show different zero-temperature phases, 
each described by a proper vector order parameter. The 
possibility of tuning a number of system’s parameters, 
in particular the interaction strength through Feshbach 
resonances, make such systems ideal to study the struc¬ 
ture of the various phases and the nature of the phase 
transitions. 

One of the easiest, but still intriguing realisation is 
represented by a two-component BEC, also known as a 
spinor condensate. Spinor condensates allow to address 
many interesting phenomena from Andreev-Bashkin ef¬ 
fect mia and fast decay of persistent currents [3], to 
(internal) Josephson effect [ll[5], or Schrodinger-cat- and 
twin-Fock-like states mz], from dimerised vortices IHl- 
to the study of quenching in classical bifurcations 
[TTI - [T8] . They represent also the basis for most of the 
recent realisation of artificial gauges in cold gases [14] . 

In this paper we specifically consider a zero- 
temperature trapped two-component BEC with an ex¬ 
ternal field that drives the population transfer (spin¬ 
flipping) between the two atomic levels (see Sec. 0 
forming the condensate. It is common to refer to the 
interconversion term as a Rabi coupling. The system 
is indeed a generalisation to non-linear atom optics of 
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the well-known linear Rabi problem and in general is 
an extension of quantum optics concepts to condensates 
[I5lll6|. It is the interplay between the intra- and inter¬ 
species two-body interaction strengths and the Rabi cou¬ 
pling strength that makes the physics of the system very 
reach. The Rabi coupling - which acts as a operators 
on each atom - tries to create an equal superposition 
of the two possible internal levels. On the other hand 
differences in the three possible atom-atom interaction 
strengths tries to favour a situation where the popula¬ 
tion of the two internal levels is unbalanced. It turns 
out that the system exhibits a second order phase tran¬ 
sition, a classical bifurcation at the mean-field level (see, 
e.g., [nillHj and in particular the experiment |S|), which 
is analogous to the mean-field ferromagnetic transition 
of the Ising model in transverse field. Moreover if the 
two component feel a different external potentials the in¬ 
ternal and external degrees of freedom are inseparable 
leading to interesting spin-orbit coupled dynamics as it 
has already been shown some years ago in mm- 

In the following we show that the static and the dy¬ 
namic response to an out-of-phase (spin) dipole pertur¬ 
bation is very rich and captures many relevant phenom¬ 
ena related to the paramagnetic-ferromagnetic-like phase 
transition of the system. A relative component perturba¬ 
tion is accessible in cold gases by applying different trap¬ 
ping potentials for different atomic internal levels. The 
spin dipole configuration is realised by applying trapping 
potentials that have the same shape, but that are dis¬ 
placed for the two components of the gas. The dynamics 
is obtained by monitoring the gas after the displacement 
is suddenly set to zero. Notice that in (TO] |20| a similar 
situation has already been realised, where instead the ex¬ 
ternal potentials were held fixed and the Rabi coupling 
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suddenly turned on. 

The main results of our analysis can be summarised as 
follow: 

(i) in the region before the bifurcation occurs, i.e., in the 
paramagnetic phase the system exhibits well defined out- 
of-phase oscillations around the equilibrium position in 
the new trapping potential. The oscillation frequency is 
in good agreement with a sum-rule approach calculation. 
The latter allows us to identify the main quantity deter¬ 
mining the spin dipole mode frequency and its relation 
with the susceptibility of the system. In the case of equal 
interaction strength the sum rule give an exact result. 
The latter is practically twice the Rabi coupling, i.e., the 
main contribution is not proportional, as it is usually the 
case, to the harmonic trapping frequency. This effect 
can be traced back in the modification of the /-sum rule, 
which is eventually due to the absence of relative number 
conservation. 


(ii) in the broken Z 2 , i.e., ferromagnetic phase the sit¬ 
uation is very different. The response of the system to 
the spin dipole perturbation is not linear and therefore 
the initial state in the displaced potentials is far from the 
equilibrium state when the potentials are the same. In 
particular the initial configuration shows a polarisation 
domain wall at the center of the cloud, but zero global 
polarisation, while in the new equilibrium state it will 
show a symmetric structure with a global polarisation. 
The dynamics is highly non-linear. After a certain pe¬ 
riod - in which the system is trapped in the domain wall 
configuration - the cloud is able to reach quickly the new 
equilibrium by spontaneously selecting one of the two 
possible polarisation. The excess energy of the initial 
configuration gives rise to small ripples in the cloud. 

It is worth mentioning here that the very same mean- 
field description we use in the following (see Sec. 0 can 
be applied, in certain regimes, to describe polariton sys¬ 
tems where the role of the polarisation is relevant (see 
e-g-, [Him]), as well as some properties of type-1.5 su¬ 
perconductors (see |23j and reference therein). 

The paper is organised as it follows. In Sec. |ll] we 
introduce the system and its description in terms of two 
coupled Gross-Pitaevskii (GP) equations. We revisit the 
emergence of a paramagnetic/ferromagnetic like transi¬ 
tion and the effect of the external harmonic trapping 
potential. In Sec. m we study the effect of a spin- 
dependent potential and the role of the spin-dipole sus¬ 
ceptibility. The latter is shown to bear a clear signature 
of the phase transition. Then, we address the problem of 
the dynamics of the spin-dipole mode both in the para- 
and in the ferromagnetic phase. In the former case (Sec. 


IV A) linear response theory combined with a sum-rule 


approach provides an accurate estimate of the spin-dipole 
mode frequency, which well compares with the numerical 
solution of the GP equation. In the ferromagnetic case 
(Sec. IV BI we show that the system exhibits ground 
state selection, after a waiting time in which the system is 
unable to leave the initial domain wall configuration. We 
found phenomenologically that this characteristic time is 


proportional to the square root of the domain wall en¬ 
ergy. 


II. GROSS-PITAEVSKII EQUATION FOR 
COHERENTLY COUPLED DECS 


We consider an atomic Bose gas at zero temperature, 
where each atom of mass m has two internal levels |a) 
and 16). The latter are typically magnetically trappable 
hyperfine levels of ®’^Rb, like |a) = |F = l,mF = — 1 ) 
(|F = l,mF = 1 )) and |5) = |F = 2,tof = 1 ) 
(|F = 2,mF = —1))- An external field is applied that 
coupled the |a) to the | 6 ) state via usually a two-photon 
transition, characterised by a Rabi splitting fl. At the 
densities of ultra-cold gases the atomic interactions are 
simply described by a contact potential with a strength 
proportional to the s-wave scattering length. For a spinor 
system three scattering lengths, Uaa, abb and Uab are 
present describing the intra- and the inter-species col¬ 
lisions, respectively. Finally the condensed phase for a 
two-component Bose gas is described by a complex spinor 
order parameter (/>a(r,t),' 0 f,(r,t)), where ' 0 i, i G {a,b} 
is the wave function macroscopically occupied by atoms 
in the internal state |z). The latter is normalised to the 
total number of atoms A/ in the state |i). The dynamics 
of the order parameter is determined by coupled Gross- 
Pitaevskii equations naEii 
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where the couplings 5 ^, with i G {a,b,ab}, are the intra- 
and interspecies atomic interaction strengths and are 
given by gi ex Oi [23] , and Va and 14 are the external trap¬ 
ping potentials. We consider the confinement to be har¬ 
monic, which is the most relevant and typical experimen¬ 
tal situation. In the following, if not differently specified, 
we consider ga = gb = g. Due to the presence of the Rabi 
coupling only the total number of atoms N = Na + Nb 
is conserved, but not its polarisation P = {Na — Nb)/N. 
The (gauge) symmetry of the system is therefore reduced 
from U{1) X 17(1) to U{1) x Z 2 , leading for an homoge¬ 
neous condensate to a gapless density or in-phase mode 
- Goldstone mode of the broken 17(1) symmetry - and 
a gapped spin or out-of phase mode (see, e.g., dZlIIH]). 
Depending on the interaction strengths and the Rabi cou¬ 
pling the ground state can also spontaneously breaks the 
Z 2 symmetry leading to P 0 0 . 

In order to describe the ground state we write as usual 
the condensate wave function as density and phase ipi = 
and use local density approximation (LDA), i.e., 
neglecting the gradient term, also known as quantum 
pressure, in Eqs. ([T). The time derivative on the LHS of 
Eqs. Q is replaced by the chemical potential whose 
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FIG. 1: Density profiles within Thomas-Fermi approximation 
for an harmonic trapping: Ua (dashed blue), nj, (light grey) 
and Ua + rif, (black) for Qablg = 1-3 and = 0.1. 


value will be fixed by requiring a total number of particle 
N. Notice that in absence of one can have two different 
chemical potentials reflecting that also P is fixed. The 
Rabi coupling originates a term of the form 17 cos{(f>a — 4’b) 
for the energy. Without any loss of generality we also as¬ 
sumed 17 to be real and positive, which fixes the phases 
in the ground state to satisfy 4>- = (pa — 4>b = Fi¬ 
nally one finds that the densities of the two component 
obey the relations (see, e.g. the review [5^ and reference 
therein) 


[g- 9ab+ , ^ ) {na - rib) = Vb- Va, (2) 

V y/rianbj 

(g + gab -^—) {ria + rib) = 2 ^ - (14 -b Va) , (3) 

V y/riaUbJ 

For the sake of simplicity and clarity, we consider a 
mean-field one-dimensional situation. The latter is ex¬ 
perimentally realised by making two of the trapping fre¬ 
quencies strong enough in order for the motion along such 
directions to be frozen. The coupling constants are in this 
case renormalised and can be simply related to the scat¬ 
tering length and the trapping transverse frequency uj± 
by gi = 2fujjj_ai for i G {a,b,ab}. It is worth noticing 
that our results do not qualitatively change in the two- 
or three-dimensional case. 

From Eq. ^ is it clear that, for equal potentials, 
Va = Vb, the system can sustain a finite polarisation only 
if gab is sufficiently large. In that case it turns out that 
the P ^ 0 states are the ground states of the system. No¬ 
tice that both the critical value of gab and P are density 
dependent. It is easy to find that the points Xp at which 
the polarised phase can exist is fixed by the condition 

gab> g + 2n/n{Xp) (4) 

with n(x) = ria(x) + nb{x) the total local density. Since 
in the harmonic trap the density decreases going outward 


from the trap center the system can exhibit two different 
regions: unpolarised tails with Ua = rib and a polarised 
core with ^ rib- Clearly, if the condition Eq. Q is not 
satisfied at the center of the trap, where the total density 
is maximum, then the whole system is unpolarised. This 
allows us to introduce a critical value of Rabi coupling 
defined by 


^cr = -jri{x = Q){gab - g) 


(5) 


For values 17 > 17cr the system is unpolarised everywhere. 
Writing 14 = 14 = rnufj^^x'^/2, the density profile Ua = rib 
is easily obtained from Eq.(|^: 


na,b{^) = no{x) 


g + n P \ 

g + gab V ) 
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where we introduced the so-called Thomas-Fermi radius 
= 2{g,+Vl)/{mujho'^) [21] and the chemical potential 
can be written as 


M = 


3 
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In the case 17 < 17cr a typical configuration within EDA 
is shown in Fig. [T] 

Let us here briefly remind that for a Bose-Bose mixture 
in the absence of Rabi coupling (17 = 0), where the rela¬ 
tive particle number can be chosen at will, the situation 
is very different. In that case there exists a first order 
phase transition to a phase separated state once gab > 9 
and the system in the trap is formed by two distinct re¬ 
gion of only one of the two component of the gas (see for 
a detailed discussion , e.g. [26ll28] ). An example of the 
structure for an equal number of atoms in both hyperfine 
levels is shown in Fig. |(al). 


III. STATIC DIPOLE POLARIZABILITY 

In this Section we calculate the static response of a 
trapped spinor gas to a spin-dipole perturbation. A 
spin-dipole perturbation corresponds to a shift of the 
harmonic traps for the two components by a quantity 
d <C Xho with Xho = \/b,/{muJho), 

Va,b = ^rnujho'^ (x ± d)^ 

= ^rnujlaX^ ± mulaxd + O(d^), (8) 

where the plus sign is for particles of component a and 
the minus one for those of component b. In the case of 
hyperfine atomic levels the displacement can be realised 
by adding a magnetic field gradient to the harmonic po¬ 
tential. 

The GP ground state solution for the spinor gas in the 
displaced potentials is reported in Fig. (see also [15]). 
where for the sake of concreteness we assume gab > 9 
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FIG. 2: Ground states of a trapped coherently-coupled Bose gas for different values of potential displacement d and Rabi 
coupling fl. We use coupling constant strength gab/g = 1.1 and g/{huihoXho) = 5. Plots (al)-(a4) correspond to djxho = 0, 
while plots (bl)-(b4) to djxuo = 0.05. For the ground state in panels 2, 3, 4 we use the values Q.IQ.cr = 0.31, 0.93, 3.11, 
respectively. For comparison we report also the ground state for a Bose-Bose mixture, i.e., fl = 0. In the latter case, being the 
number of particle in each component fixed, no global polarisation appears, and the ground state (al) and (bl) are essentially 
equal. The effect of quantum pressure can be clearly noticed in the plots (a2) and (a3) (analogues to Fig. [^. The bifurcation 
points are not sharp as instead predicted by the Thomas-Fermi approximation. 


to show the difference between a mixture and a coher¬ 
ently driven spinor gas. For d = 0 (row a) we see the 
features of the f2-induced phase transition: below the 
critical value the linear coupling prevents the phase sep¬ 
aration by creating a global polarisation in the system 
(Fig. [^plots (a2) and (a3)), while a mixture without any 
Rabi coupling = 0 would be in a phase separated state 
(Fig. I plots (al)). Above the critical value the gas is 
unpolarised (Fig. plot (a4)). Applying a potential dis¬ 
placement (row bjmakes the local polarisation different 
from zero as shown in Eq. In this case even a small 
potential difference makes the ferromagnetic part of the 
gas strongly polarised and as a result a magnetic domain 
wall is created at the center of the trap. 

In order to calculate the spin-dipole susceptibility we 
first determine the spin-dipole moment D, defined as 

D=^J x[na{x) - nb{x)]dx. (9) 

The spin-dipole susceptibility is then defined by the limit 
Xsd = lim D/X (10) 


where A = dmujf^^ is the perturbation associated with the 
spin-dependent component of the potential (|^. 

In the global paramagnetic phase (17 > flcr) it is easy 
to obtain an analytical expression for Xsd within LDA. 
In this case one can employ the energy functional 

[x'^^ina-mf - Xx{na-nt)] dx (11) 


relative to the spin degrees of freedom of the problem, 
where 


Xs = - . n/ 

g- 9ab + il/no 

is the spin (magnetic) susceptibility for an homogeneous 
system of density 2no (see, e.g., [IS]). Variation with 
respect to the spin density (ua — rib) yields the result 


na{x) - nb{x) = xXxs{no{x)), (13) 


and the spin dipole polarizability finally reads 


Xsd 



(14) 


After integratiion we obtain the result 

D _ 9 T gab j. ( __ 

A 9-9ab[ \{g-9ab)no 

for the dimensionless ratio D/d = moj/^^Xsd where 
we have introduced the dimensional function f{a) = 
3a (l — ^/l + a arccoth(vT+^)) [STj and used the no¬ 
tation no = no(0). 

A couple of comments are due here. First of all let us 
consider the case of a Bose-Bose mixture, i.e., 17 —t 0. In 
this case /(a —?► 0) —t 0 and the spin-dipole susceptibil¬ 
ity is simply proportional to the magnetic susceptibility 
Eq. ( |I2| ). Therefore also Xsd diverge at the (miscible- 
immiscible) transition point gab —t g~ [38]. Physically 
this is due to the fact the two gases become globally im¬ 
miscible at the transition point since the latter condition 
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FIG. 3: Dipole D as a function of traps displacement d for 
different values of D/Dcr and for gabjg = 1 .1 a s in Fig. 
Dotted lines are analytical results from eq. (151 for the four 
bigger values of fl, points are numerical data and grey full 
lines are only a guide for the eyes. 


is density independent. As we will see in the next Sec¬ 
tion the divergence of Xsd leads to a zero frequency (soft) 
spin-dipole mode. 

On the other hand for finite O the paramagnetic- 
ferromagnetic transition point, namely gab = 9 + ^/no, 
depends on the density. Therefore the spinor gas starts 
becoming ferromagnetic at the center of the trap only. 
The quantity Xsd, being density integrated, remains fi¬ 
nite at the transition point (indeed /(—I) = —3) leading 
(see next Section) to a finite frequency for the spin-dipole 
mode. This behaviour is very general and it has already 
been pointed out for the Stoner (or itinerant ferromag¬ 
netic) instability in the context of cold gases by two of 
us [2S]. 

Above the critical point the response of the system is 
no longer linear. The system is partially ferromagnetic 
and has the tendency to form a magnetic domain wall at 
the centre of the trap (see Appendix). 

A detailed analysis of the behaviour of D/d is shown 
in Fig. where we calculate numerically the spin-dipole 
of the gas as a function of the trap separation d with the 
choice gab/9 = 1.1. Above the critical Rabi frequency we 
see that indeed linear response applies and the analytical 
expression Eq. (15) works very well. Notice that the 
spin-dipole moment allows for a clear identification of the 
phase transition point, above which the induced dipole 
moment D changes its behavior as a function of d. 


IV. SPIN DIPOLE DYNAMICS 


In this section we study the dynamics of the system. In 
particular we prepare the system initially in the ground 
state of very slightly displaced external potential and 


then suddenly set the displacement to zero. As one can 
expect the physics is completely different depending on 
whether the system is paramagnetic or ferromagnetic. In 
the earlier case the system shows a well defined out-of- 
phase oscillation, the spin-dipole mode. The previously 
calculated spin-dipole polarizability plays an important 
role in characterising the behavior of the spin-dipole fre¬ 
quency [ini 1301 • Notice in particular that for two inde¬ 
pendent condensates (fl = gab = 0) the spin-dipole fre¬ 
quency would simply coincides with the trap frequency 
ujho- In the ferromagnetic case the system evolves accord¬ 
ing to a highly non-linear dynamics and it shows ground 
state selection. We analyse the two cases separately in 
the next two sections. Some details on the numerical 
solution of the GP equations can be found in Appendix 
B and reference therein. For the interested reader we 
include in the Supplementary Material the real time evo¬ 
lution of the system in different regimes. 


A. Paramagnetic phase: sum rule approach 

In the paramagnetic phase, as shown in Fig. § (a4)- 
(b4), a small trap displacement corresponds to a small 
deviation with respect to the ground state at zero dis¬ 
placement and therefore linear response theory can be 
applied. The dynamics we consider in this case coin¬ 
cides with the dynamical response of the spinor gas to 
the spin-dipole operator Sd = A very pow¬ 

erful tool to estimate the frequency of collective modes is 
the so-called sum rule approach [STJ [32] . This approach 
has been very successfully employed for the dynamics of 
both cold gases and nuclei. We remind here simply that 
sum rules are defined for an operator F as 

mk{F) = Y,mF\n)\^iE^-Eo)^ (16) 

n 

and they represents the moments of the strength distri¬ 
bution function relative to F. The sum rule approach has 
the merit of providing a direct way to obtain an upper 
bound of collective mode frequency through the ratio of 
different sum rules, and therefore giving an understand¬ 
ing of the collective mode frequency in terms of static 
macroscopic quantities [32) . 

In our case the operator of interest is Sd and we we 
use the energy weighted and inverse energy weighted sum 
rule, i.e., 




mijSd) 

m-i{Sd)' 


(17) 


They are particularly suited in our case. The energy 
weighted one (mi) is easily rewritten in terms of a double 
commutator as mi = (I/2)(0|[S'd, [i7, S'd]]|0). The only 
terms in Ff that do not commute with Sd are the kinetic 
energy and the Rabi coupling Ff^ = The 

former gives the usual Nti^ / (/2m) contribution, while the 
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FIG. 4: Spin-dipole frequency as a function of Q for differ¬ 
ent values of interactions. Lines are analytical results from 
Eq. ( |20[ ) and points are numerical data. In order to have a 
fully paramagnetic phase for gab > g one needs fl > flar (see 
text and Eq. ([^), the value of flcr for gabjg = 1.02 (blue di¬ 
amonds) is shown by the dashed line. In the Supplementary 
Material, two videos show the oscillations of the clouds for 
dab!g ~ 0.9 and gabjg ~ 1 (both at If = O.Bhujho)'- 1.0 sec of 
the video corresponds to l.Ocuhot. 



gab/g 


FIG. 5: Spin-dipole frequency for a Bose-Bose mixture, i.e. 
II = 0, as a function of the ratio gab/g- Line is the analyt¬ 
ical result and points are numerical data. In this case the 
spin-dipole frequency goes to zero at the phase separation 
transition point. 


law 

LUsoi^ = 0 ) = LOho\l ^ , ( 21 ) 

V 9 + dab 


latter is straightforwardly evaluated as —AVlx^ax- Aver¬ 
aging on the ground state, we obtain the result 

/*^TF 

mi = N -1-811 / x'^nQ{x)dx. (18) 

2 m Jq 

The inverse energy weighted sum rule (m_i) is directly 
related to the susceptibility of the ground state through 
the relation 


N 

m-i = —Xsd, 


(19) 


and using the definition Eq. (10) together with the result 


Eq. (141 we obtain the following upper bound to the 


spin-dipole frequency 


2 _ 2f 9-9ab \ l ^ + 8flno{g + gab)/i5hujiy 

WSD OJho + [ 1 + 

( 20 ) 

Notice that the equality in Eq. 0 is attained when 
the whole strength is exhausted by a single state. 

In Fig. 1^ the sum-rule result is compared with the 
predictions of the solutions of a time dependent Gross- 
Pitaevskii calculation. As already mentioned from the 
numerical or experimental point of view, the excitation 
of the spin-dipole mode is achieved starting with an equi¬ 
librium configuration in the presence of slightly displaced 
trapping potentials, as described by Eq. (|^, and sud¬ 
denly setting d = 0. 

Notice that at the transition point the frequency does 
not go to zero, since for the reasons explained in the pre¬ 
vious Section Xsd (or m_i) does not diverge at that point. 
This has to be compared with the mixture case, which is 
recovered sending 12 —>■ 0. In this case the spin-dipole fre¬ 
quency vanishes close to the critical point following the 


and the sum-rule approach give the exact result as shown 
in Fig. 

Sum-rules give the exact result also for the intrinsic 
SU{2) symmetric point gab = 9 (and 12 yf 0 in gen- 
eral)and gab = g as it can be seen in Fig. (red tri¬ 
angles). The magnetic energy of the spinor gas in this 
regime depends on the relative density only through the 
Rabi coupling which breaks the SU{2) symmetry of the 
system. The spin-dipole frequency behaves in this case 
as 




( 22 ) 


which is essentially twice the Rabi frequency and there¬ 
fore almost independent of the tapping frequency. The 
latter unusual result for a trapped gas is due to the corre¬ 
lation between the internal and external degrees of free¬ 
dom that in particular lead to the modification of the 


/-sum rule, see Eq. (18). 


In the more general case, when both 12 and (gab—g) are 
different from zero, the frequency is given by the full eq. 
(20) in which both the coherent and the interspecies s- 


wave couplings play a role. In this more general case one 
observes that that the sum rule approach provides only 
an upper bound to the numerical solution, due to the 
appearance of more frequencies in the numerical signal 
resulting in beating effects. 

In the Supplementary Material we include two videos 
showing the oscillations of the clouds in the paramagnetic 
case for gab/g = 0.9 and gab/g = 1 and 12 = O.'bhiOho- The 
real time evolution shows clearly the presence of only one 
frequency in the intrinsic SU (2) symmetric case and the 
appearance of more frequencies when the gab 7 ^ 9- 
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FIG. 6: (a) and (b): real time spin-dipole oscillation (black solid lines), polarization of the system (blue solid lines) and ground 
state polarization values (dashed lines); parameters are gab = l-lfl, g/ifvjJhoXho) = 5 and Q.lhiiJho = 1-5 (a), fl/hcuho = 2 (b). 
In the Supplementary Material a video shows the oscillation of the clonds of panel (a), 1.0 sec of the video corresponds to 
l.QuJhot. (c): waiting time as a function of 11, points are numerical data and line is a fit of data with function A^/a where a is 
the surface tension of Eq. (231. 


B. Ferromagnetic phase: ground state relaxation 


In the previous section we have studied the dynamics 
for a completely paramagnetic gas, i.e., 17 > Her- The 
behaviour is very different when the system presents a 
ferromagnetic behavior. In this case the ground state of 
the system with equal trapping potentials is polarized as 
shown in Fig. [^plots (a2) and (a3). When the traps are 
shifted, the ground state is instead globally unpolarized 
{Na = Nb) but with a large spin-dipole moment (depend¬ 
ing on the values of Q and d) as one can see in Fig.j^plots 
(b2) and (b3). Therefore the initial state and the ground 
state are very far from each other. This circumstance 
results in a non-trivial non-linear dynamics as shown by 
the dynamics of the spin-dipole and of the polarisation 
reported in Fig. At the beginning, the spinor gas 
oscillates around the initial configuration, trapped in the 
unpolarized state. After a certain time, Twait, the domain 
wall starts moving and a finite polarisation appears. The 
system then bounces back and forth between the initial 
magnetic state and its magnetic ground state to even¬ 
tually relax to the latter one. [39] An example of such 
dynamics can be viewed in the video in the Supplemen¬ 
tary Material. If the global polarisation of the ground 
state is large, the effects of non linearity and the number 
of bounces are large. When the system is slightly in the 
ferromagnetic regime no bounces are observed and the 
system after Twait soon reaches its ground state (see right 
lower panel in Fig. |^. Notice that even if the system is 
isolated, it can approach in the long time limit an asymp¬ 
totic steady state, as a result of destructive interference 
of several time oscillating factors, present in the evolution 


of expectation values of observables. Specifically, in the 
case of a large and dense collection of frequencies, the in¬ 
terference phenomenon results in a dephasing mechanism 
similar to inhomogeneous dephasing. 

The total energy of the system is still conserved, the 
algorithm used (see Appendix]^ does not contains any 
dissipative mechanism and we explicitly check that the 
total energy does not change during the evolution. At the 
end of the real time evolution we get the ground state pro¬ 
file with superposed some high frequency perturbations 
carrying the extra energy. 

As we already mentioned in Sec. II the initial configu¬ 
ration in the ferromagnetic case contains a domain wall 
at the centre of the trap. We have identified a close re¬ 
lation between the observed waiting time and the square 
root of the domain wall energy (see Appendix [A|) 


l(5-5a&)« + 2np/2 


(23) 


From an intuitive point of view, the higher is the energy 
of the domain wall, a, the more time is required to the 
system to relax from the kink into one of the ground 
states of the system; accordingly, it is expected a rela¬ 
tion of proportionality between the waiting time and a. 
A standard field theoretical estimate of the average tun¬ 
nelling time cannot be straightforwardly performed since 
only close the transition our field theory resembles an or¬ 
dinary (/)^ theory (see Appendix A for details); for this 
reason, we benefited of a numerical fit to extract with 
surprising accuracy the relation, Twait oc ^/a^, as shown 
in Fig. [^ 

The fact that for 17 —>■ 0 the waiting time diverges can 
be easily understood noticing that the initial state and 















the ground state are very far from each other (see, e.g., 
panels (a2) and (b2) of Fig. [^. Eventually, in the strict 
n = 0 case, the system cannot reach the totally polarised 
ground state and it remains in the phase separated state 
(see panels (al) and (bl) of Fig. |^. 


V. CONCLUSIONS 

In the present work we analyse in details the static 
and dynamic response of a trapped coherently driven 2- 
component condensate to spin-dipole probe. We show 
that the spin-dipole susceptibility is a good quantity able 
to identify the appearance of a ferromagnetic-like region 
in the cloud. 

For the dynamics we study the spin-dipole mode fre¬ 
quency by starting in a configuration with displaced har¬ 
monic potentials, which are suddenly brought to the same 
value. When the system is paramagnetic such a fre¬ 
quency is well reproduced by a sum-rule approach. In 
particular the /-sum rule is strongly modified by the Rabi 
coupling in the symmetric interaction case {g = gab) and 
the inverse energy sum rule is proportional to the spin- 
dipole susceptibility and coincides with the second spa¬ 
tial momentum of the local magnetic susceptibility (see 
Eq.@). 

When the system has a ferromagnetic domain a linear 
response cannot be applied anymore and the dynamics 
is highly non-linear. The initial configuration within dis¬ 
played potentials is unpolarised and contains a magnetic 
kink centred at the origin. The dynamics is trapped for 
a time, Twait, in the initial configuration, after which the 
system is able to relax to its polarised ground state. We 
find that Twait is proportional to the square root of the 
kink surface tension. 

Our study improves the characterisation of coherently 
driven BECs, enlightening their differences with respect 
to Bose-Bose mixtures. Moreover measuring the spin- 
dipole dynamics opens new perspective to experimentally 
access important magnetic properties of the system, as, 
e.g., its susceptibility or the domain wall surface tension. 


where the first term arises from the kinetic energy and 
the term 


W{M) 


4 


[g - 9ab)M‘^ - fInVl - M2, 


(A2) 


accounts for the density-density interaction and the Rabi 
terms. For a homogeneous magnetization minimisation 
6E/6M = 0 leads to the usual equation for the para- 
and ferro-magnetic like states. From Eq. (Al) one sees 


that close to the phase transition, i.e., M 1 a. standard 
Ginzburg-Landau theory for the order parameter M, is 
valid, where the kinetic energy is just the the square of 
the gradient of M and the effective potential takes the 
usual quadratic plus quartic form 


W{M) = 


g- gab + 


2Vl 




nU, 




= 

2 




(A3) 


As usual the Z 2 symmetry broken ground state is ob¬ 
tained for r < 0. A kink in M is the field solution 
interpolating between the two degenerate minima. Its 
surface tension, a, which coincides with its energy in a 
one-dimensional situation, can be easily computed |34] 
yielding the result 


cr oc 


|j,|3/2 / h^n'^ \Sgn -|- 20^/^ 


oc 


m 


m 


n 


(A4) 


Appendix B: Numerical method 

All numerical data presented in this paper have been 
obtained solving the GP coupled equations by means of 
the split-operator method and by treating the kinetic 
term in Fourier space [351 I36j . The initial wave func¬ 
tions 'ipaix,t) and are evolved for a time step At 

alternately by the kinetic, potential and Rabi terms of 
Hamiltonians in Eq. ([^: 


Appendix A: Magnetic domain wall surface tension 


In this Appendix we briefly show how to approxi¬ 
mate the energy functional for the magnetisation for a 
spinor condensate with a classical one dimensional (jA 
(Ginzburg-Landau for the phase transition) theory jddj . 
From the latter it is easier to show the existence of a 
kink or domain wall in the magnetization, and we com¬ 
pute in this regime its surface tension. In the symmetric 
case ga = gb = g and considering an uniform total den¬ 
sity n = Ua + rib, the relative density or magnetisation 
M = {ria — nh)/n, enters in the energy density 


E{M) 



tAn{VM)^ 
8to(1 - M2) 


+ W{M) , 


(Al) 


ipi{x,t) ^ cosh{—nAt)Ipi{x,t) + sm'h{—ilAt)ipj{x,t) 


where i = a,b ^ j and Eq. (Bl| is for the imaginary 
time evolution. One can simply obtain the same set of 
equations for real time evolution by changing At in iAt. 

An algorithm of this type is symplectic, this means 
that the method exactly simulates a Hamiltonian 
with Hai — Ha power series in At. Advantages of using 
symplectic integrators are that there is no drift in energy 
due to the exact conservation of Hai and that the phase- 
space volume is exactly conserved. 

In order to obtain the ground states we ran the imag¬ 
inary time evolution starting from an initial trial wave 
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function built from both random density and phase dis¬ 
tributions, in order to prevent the algorithm from reach¬ 
ing false metastable states. For the dynamics we loaded 
the ground states obtained with displaced traps and let 
them evolve using the same algorithm but with real time 
and with the equal trapping potential Hamiltonian. The 
values of polarization P and of spin-dipole moment D are 
calculated and saved at each time step and then analysed 
to obtain the frequencies. This last step is not always 
so straightforward, sometimes the signal contains more 
than one frequency and damping occurs. In this cases 
we perform a Fourier analysis of the data and we keep 
the maximum-amplitude frequency. 
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